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Abstract. A homogeneous set of UBV photometry (354 
data points obtained between 1983 and 1998) for the 
Be/X-ray binary A0535+26 = V725 Tau is analysed, aim¬ 
ing to look for possible periodic component (s). After sub¬ 
traction of the long-term variation it was found that only 
a ft 103 d periodic component remains in the power spec¬ 
tra in both the V and B colour bands. The probability 
of chance occurrence of such a peak is less than 0.1%. 
There are no signs of optical variability at the X-ray pe¬ 
riod (« III d ). We discuss possible reasons for a 103-day 
modulation and suggest that it corresponds to a beat fre¬ 
quency of the orbital period of the neutron star and the 
precession period (ft 1400 d ) either of an accretion disc 
around the neutron star or a warped decretion disc around 
the Be star. 

Key words: stars: binaries: general - stars: circumstellar 
matter - stars: individual: V725 Tau - stars: pulsars: in¬ 
dividual: A0535+26 - stars: variables: general - accretion, 
accretion discs 
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1. Introduction 

The X-ray binary A0535+26 has been an object of in¬ 
terest to both observers and theorists for a quarter of a 
century, from the moment of its first documented X-ray 
outburst in 1975 (Rosenberg et al. 1975). As soon as an 
optical counterpart was identified almost simultaneously 
by several authors (Liller 1975| ; Murdin 1975| ) with the 9th 
magnitude Be star HDE 245770, later named V725 Tau, 
there were several attempts to find the orbital period of 
the system, using both X-ray and optical (spectral and 
photometric) data. 


Priedhorsky & Terrell (1983) have shown that all but 
the most powerful X-ray outbursts (1975 April and 1980 
October) take place with a period «111 days. Now the 
general opinion is that the orbital period is 110-1 ll d . The 
most reliable determination of the X-ray period and initial 
epoch is currently that of Motch et al. ( 1991 ): Pxr = 
lir?38 ± 0^11 and T 0 = JD2446734.3 ± 2‘?6. 

This and other periodicities were claimed to have been 
found in optical (photometric and spectral) data. For in¬ 
stance, Hutchings et al. ( 1978| ) found that the radial ve¬ 
locities of absorption lines of HDE 245770 are variable, 
and suggested several probable periods: 28.6, 48 and 94 
days. Some years later and based on a very similar dataset, 
Hutch ings (1984) found a 112 day period. Guarnieri et al. 
( 1982 ) noted at first possible 32, 63 or 77 day peri ods in 
their photometric data, but later (Guarnieri et al. 1985 ) 
found modulation in the V band photometry with the 
proposed orbital period of 110 days. Gnedin et al. ( 1988 ) 
made a Fourier analysis of V band photometry obtained 
during 1981-1985 and marked out periods of 1100, 103 
and 28 days, but found no periodicity corresponding to 
the X-ray period. 

Besides the determinations noted above, there are sev¬ 
eral more publicati ons o n this subject, listed by Giovan- 
nelli and Graziati (1992). 


830 and 507 days. However their light curve modeling does 
not agree with the observational data and their prognosis 
of the photometric behavior of HDE 245570 for 1994- 1997 
was not confirmed by later observations (see Fig. 1 in Hao 
et al. 199C| . Fig. 1 in Lyuty & Zaitseva 2000| and Fi g. p| in 
this work). It should be noted also that in Hao et al. (1996) 
no correction has been made for different zero-points and 
colour systems of the data sampled by different groups. 

As soon as our data set had grown su bstant ially 
and surpassed that used in Gnedin et al. ( 1988 ), as 
well as other published data sets, we attempted to 
search for periodic components in the light curves of 
A0535+26/V725 Tau. In Sect. 2 we describe our tech¬ 
niques of data analysis and obtain a value of the optical 
period, distinctly different from the orbital one, while in 
Sect. 3 we discuss possible reasons of the light modulation 
and suggest that it is caused by interplay of the orbital and 
precession motions. In Sect. 4 we summarize the results. 


2. An analysis of optical data 

All the photometric data discussed herein have been ob¬ 
tained in the period 1983-1998 with a photon-counting 
photometer at the 60 cm telescope of the GAISh Crimean 


laboratory (for details see Lyuty & Zaitseva 2000). More 
than 350 individual UBV measurements have been made; 
the photometric accuracy is 0“005 in B, V and 01 n 008 in 
U. As the long-term behavior of V725 Tau has been dis¬ 


cussed inseparate papers (Clark et al. 199E, Lyuty & Za¬ 
itseva 200C ), we restrict ourselves here only to the analysis 
of variability on the time scales close to the X-ray period. 

We calculate the power spectrum and spectral window 
of the data set as 


In a recent paper by Hao et al. ( 1996 ) all the published 
data sets from 1981 to 1993 have been analysed together. 
Their analysis have shown only two significant periods - 
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JD 2400000 + 

Fig. 2. a The light curve of V725 Tau in V band, b 
smoothed light curve with a 50 d window, and c residu¬ 
als after subtraction of the smoothed curve. 


Fig. 1. a Power spectrum and b spectral window in the 
V band for our 1983-1998 observations after subtraction 
of the linear trend. The highest peak corresponds to a 
(quasi)period of ~ 1400 d . 


respectively, where i = V"-I, tj denotes the time of ob¬ 
servation, and uij and to, correspondingly, the individual 
and mean values of the brightness. 

Fig. 0a shows the power spectrum of the V band data 
(linear trend removed) for the period from 1983 to 1998. 
The largest peak in the spectrum corresponds to ~ 1400 
days and is readily explainable: a wave with that charac¬ 
teristic time scale is easily seen in Fig. 0, superimposed 
on a smooth decay in brightness, in the optical and in¬ 
frared light curves (see also Larionov 1993| and Clark et 
al. 1999). In the region of the orbital period (~ 0.01 d -1 ), 
only a faint peak corresponding to P = 103 d is seen. At 
the frequency corresponding to a 111 day period, there is 
no peak at all. 

In order to judge with sufficient confidence the reality 
of the orbital plus close and/or related periodicities, one 
should remove the longer-period component (s) of variabil¬ 
ity. In our case it is not sufficient to subtract the linear 
trend, and we need to take into account the slow light 
variations with a (quasi)period of ~ 1400 days. We have 
constructed an approximating data set using the method 
of a sliding mean with the window value A, replacing raw 
data rrii for each time t,; by the weighted mean: 


m 


; =-2.5 log 


k 




Pi 


10 


—0.4m,- 


( 3 ) 


where k is the number of data points within interval [U — 
A, ti + A], and the weight of jth point is determined as 
Pi = exp [-(c5fj/A) 2 ] , (4) 


where Stj is the time span from the jth point to the cen¬ 
ter of the window. The optimal value of the smoothing 
interval was searched by trial and error within a range 
20 d < A < 100 d . The criterion for A selection was the 
signal-to-noise ratio of the power-spectrum peaks in the 
region of interest, i.e. around the orbital frequency. The 
raw and smoothed light curves, and also the residuals be¬ 
tween them, are plotted in Fig. 0. The power spectrum 
for the residuals rrij — to' is shown in Fig. 0. We have 
found that, independent of the value of the smoothing 
interval adopted, the most prominent peak in the power 
spectrum corresponds to a period of 102/83 ± Orl, while 
the low-frequency components are most effectively sup¬ 
pressed and the highest signal-to-noise ratio is obtained 
when A = 50 d . We have obtained an analogous result for 
the B band as well. 

The amplitude of the sine-wave obtained is OVffilS in 
V and O^Ol in B, while in U the 103-day peak is not 
seen above the noise. The significance of the detection of 
the faint sine-wave superimposed on the large-amplitude 
slow variability can be tested by adding an artificial low- 
amplitude harmonic to the raw data. We added a sine- 
wave with P = 76 d and amplitude 0W015 to the raw V 
data, then repeated the same cleaning routine as described 
above (Eqs. 0 and 0 ). Fig. 0 demonstrates that the low- 
frequency filtration method described above can be used 
to detect, at least for the data set discussed , periodic vari¬ 
ability on a time scale of ~ 100 d with amplitude > 0“01. 

In order to test whether the 103-day period is real or 
an artifact intrinsic to some part of the data set, we have 
split the initial data set into two parts corresponding to 
values above and below the linear trend and have repeated 
the same procedure. We have found that the 103-day har¬ 
monic is intrinsic to both parts of the initial dataset, al¬ 
though for the low-brightness part the significance of the 
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Fig. 3. a Power spectrum of the residuals V — V' on our 
1983-1998 data. The most prominent peak corresponds 
to P = 102/92. b To check the efficiency of smoothing, 
an artificial sine-wave with P = 76 d amplitude (F'OIS was 
added. Additional peaks in the spectrum, the most promi¬ 
nent of which are at P = 63 d , are caused by aliasing of 
the P = 76 d period with the peaks in the spectral window. 
The slight decrease in the power of P = 102/83 peak, as 
compared to the upper panel, is caused by a redistribu¬ 
tion of power as a result of the cleaning process on our 
non-uniformly spaced data set. c Power spectrum after 
subtraction of 102/83 and 91/07 period sine-waves from 
the initial light curve. 


103 d peak is slightly less. Fig. g gives the phase curves 
of residuals V — V' for the whole data set and for the 
“upper” and “lower” parts of it. It is easily seen that 
both the shape and initial phases are the same. We calcu¬ 
late the ephemeris of the small-scale optical variations as 
JD M in = 2448002(±3 d 0) + 102/83(±0/l5) • E, where E is 
the epoch number. 

The probability of chance occurrence of a peak with 
amplitude P max in the power spectrum with a mean 
value P m ean can be estimated as $ = 100% • {1 — [1 — 
exp[— (P max /f’mean)]] JVind }; for non-uniformly spaced data 
the number of independent frequencies Wnd = —6.362 + 
1.193 • N + 0.00098 • N 2 , where N is the number of obser- 



0.0 0.5 1.0 1.5 2.0 


Phase of 102 cl .83 period 

Fig. 4. a V — V' dependence on the phase of the 103- 
day period for the whole set of data, b for the “upper” 
and c for the “lower” part of it. A 0“015 sine-wave is 
superimposed on each panel. 


vations (Horne and Baliunas 1986). In our case N = 354 
and A'jnd = 539; after subtraction of the slow component 
as described above (Eqs. 1 andg), < 0.1%. 

Nevertheless, there are some peaks in the frequency 
range 0.007 — 0.017 with power > 10 -5 . In order to test 
their reality, we hav e followed the recipe given by Horne 
and Baliunas ( 1986 ) for irregularly spaced data and sub¬ 
tracted a 103 d period sine-wave from the initial light 
curve. After this we find that only a peak at a frequency 
corresponding to 91/07 and its aliases remain in the power 
spectrum. This feature replaces the 103 d periodicity from 
JD 2449400 (1993 April) onwards. It is remarkable that 
this change of periods coincides with the time of the most 
prominent X-ray outburst over the past three decades. 

Fig. 3c displays power spectrum after subtraction of a 
sine-wave of period 103 d before 1993 April and 91 d after 
that time. Comparing Fig. 3a and Fig. 3c, we conclude 
that only the 103 d and 91 d peaks are real, although the 
latter feature is present in only 1/5 of the total data set, 
thus preventing us from making exact estimates of the 
ephemeris. A slight enhancement of the noise level around 
0.01 d _1 on Fig. 3c is most probably caused by the devi¬ 
ation of the real signal from a true sine-wave. 
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3. Discussion 

We stress once again that the 111-day orbital periodic¬ 
ity is not revealed in the optical photometry. It is then 
natural to suppose that one of the constituents of the 
total radiation of the system, besides that of the op¬ 
tical star, is that of a precessing disc - either an ac¬ 
cretion disc around the neutron star or a tilted/warped 
equatorial envelope around the optical star. The period 
of precession, whichever precessing body, would then be 
Pprec = l/(l/Popt - 1/ -Px—ray) = 1360 d ± 30d, and ~ 900 d 
after 1993 April. In the following we consider these two 
scenarios separately. 


3.1. Accretion disc 


In the case of the neutron star (NS) accretion disc, the axis 
of the disc is inclined to the orbital plane of the system 
and, additionally, is counter-precessing with a ss 1360 d 
period. In this case the disc cross-section, as seen from the 
optical companion, would be changing with a frequency 
equal to the sum of orbital and precession frequencies, 
which would lead to the 103-day modulation observed. 

The existence of an accretion disc in this binary sys¬ 
tem at least during some X-ray outbursts is firmly estab¬ 
lished - it is confirmed by quasi-periodic os cillat i ons o f 
X-rays during X-ray outbursts (Finger et al. 199(j . 199G) 
and by neutron star spin-up episodes (Nagase et al. 1984 ). 
However, the question as to whether the disc exists per¬ 
manently or is formed just before the outburst remains 
unresolved. 

A similar model was proposed to describe the opti¬ 
cal light curves of some other X-ray binaries: Cen X-3, 
LMC X-4 and others (see, e.g., Heemskerk & van Paradijs 
|1989| ) - however, unlike A0535+26, these systems show or¬ 
bital variability besides precession. This can be explained 
by the fact that A0535+26 is a wide binary, and the ef¬ 
fects caused by the ellipsoidal shape of the optical com¬ 
ponent and/or its X-ray heating are unobservable. Also, 
in the case of A0535+26 the large eccentricity of the or- 
0.5 (Finger et al. 1994, 1996) - plays a major 


bit - 

role. As a result, both the effective cross-section of the 
disc to the stellar wind and its illumination in periastron 
and apastron would differ by ~ 9 times; when ss 1360 day 
precession is added, this enables the 103-days modulation 
to be observed. 

Meanwhile, the disc’s input to the total optical radia¬ 
tion of the system can be as small as a few per cent (Lyuty 
& Zaitseva |2000 ), and it would be quite a difficult task to 
distinguish it based on the difference between its photo¬ 
metric and/or spectral parameters and those of a giant 
star and its envelope. 

The disc is not expected to rotate as a solid body, and 
its structure is dependent on X-ray outbursts occurring at 
the surface of the neutron star and subsequent re-radiation 
of stored energy. Because of this, neither the amplitude of 


precession modulation nor the period and phase of preces¬ 
sion should be stable. 

Meanwhile, we see that during the first 10 years of 
our observations (up to 1993), the 103-day periodicity of 
the optical variations is retained, which means that the 
precession period is rather stable, changing from ~ 1360 d 
to ~ 900 d after the major X-ray outburst. Nevertheless, 
as expected, in the total radiation of the system the rel¬ 
ative contribution of the component connected with the 
disc precession is not constant. Fig. || illustrates the phase 
curves in V band after subtraction of the slow component, 
separately for the ascending and descending branches of 
the initial light curve shown in Fig. |[ It is evident that the 
103-day modulation is practically absent when the light 
level is increasing, and, on the contrary, is quite distinct 
when the light fades. 


-0.05 
? 0.00 
0.05 

-0.05 
? 0.00 
0.05 


Fig. 5. V — V' dependence on the phase of the 103-day 
period a for descending and b ascending branches of the 
slow component of the light curve (see Fig. |^). 
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One of the possible mechanisms of stimulation of the 
accretion disc precession is free precession of the neutron 
star. It was shown by Schwarzenberg-Czerny ( |1992 ) that 
-Ppr/-P sp in ~ 10 6 , where P pr is the period of free preces¬ 
sion of a neutron star and P sp i n is the period of its axial 


rotation. In the case of A0535+26 P st 


104 s and cor¬ 


respondingly P pr ft: 1200. It is natural to suppose that 
the proximity of this period to the characteristic time of 
the large-scale photometric variability is not coincidental. 
What then can serve as a “transmission link” from the 
freely precessing neutron star and its accretion disc to the 
equatorial envelope of the optical component? Let us be 
reminded that the neutron star in A0535+26 system has 
a powerful magnetic field, « 10 13 G. Periodic changes in 
orientation of the magnetic field relative to the equatorial 
envelope would lead to periodic changes in the strength 
of their interaction, first of all near periastron; the shock 
wave arising in the envelope would lead to (quasi)periodic 
ejections of the envelope. Clark et al. (1999) and Lyuty & 
Zaitseva (2000) argue that the long time scale variability 
of A0535+26 in the period we are analysing is explained 
by the successive expulsions of the Be star envelope that 
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occur on a characteristic time-scale of ss 1400 days. We 
have to note, however, that the detailed description of the 
interaction of the Be star envelope and NS magnetic held 
needs additional modelling work that is beyond the scope 
of this paper. 


3.2. Be-star equatorial envelope 


A quite different approach may be considered on the basis 
of recent analysis of spectral variability of another Be/X- 
ray binary, V 635 C as=4U0115+63, made b y Neg ueruela 
and Okazaki ( [2001 ) and Negueruela et al. (2001). These 
authors argue that in V635 Cas, just as in similar Be/X- 
ray systems, a major role is played by an equatorial de- 
cretion disc around the optical star. This disc is truncated 
as a result of tidal/resonant interaction with the neutron 
star companion. At a certain stage the disc becomes un¬ 
stable, tilts and warps and starts to precess. Later on, 
the disc is disrupted due to interaction with the orbiting 
neutron star, and a giant X-ray outburst may occur. The 
above authors argue that this model does not imply any 
substantial change in the optical star’s mass loss rate. 

Following that idea, we may suppose that the process¬ 
ing body in the A0535+26 system is a tilted/warped disc 
around Be star. Then it is natural that the gravitational 
pull of the neutron star in the moments of closest ap¬ 
proach of the NS to the disc causes its distortion, which is 
observed in the light curve as minor variations superposed 
on global changes. Due to precession, the times of closest 
approach do not coincide with the periastron passages, 
but rather precede them by ss 8 d after each orbit. Within 
this approach we can find a natural explanation for the 
fact that during the ascending parts of the global light 
curve, no 103-day modulation is seen: it only means that 
the decretion disc is not large enough to become tilted and 
warped and therefore it does not precess. 

To further develop this model, one may speculate that 
the changes of projection of the warped disc to the plane 
of sky, as seen from the optical companion, should cause 
substantial variations with the processional period, super¬ 
posed on the optical light curve. This is actually observed: 
we noted before the coincidence of the 1400 d time-scale of 
global optical variations with almost the same period of 
precession - whichever the precessing body. 

One might expect that the X-ray activity within the 
system may also in some way be connected with the phases 
of the 103 d period. We tried to check whether such a cor¬ 
relation exists, and found that out of 18 documented out¬ 
bursts with amplitude > 0.2 Crab (Giovannelli & Graziati, 


1992 and Finger et al. 1996), 7 peaked in a narrow interval 
0.6 — 0.7 of the optical phase curve (Figs. [| and |5), where 
phase 0.0 refers to optical minimum (see Table 3 ). This 
means that X-ray outbursts “prefer” to happen 10 to 20 
days after the optical maxima. If we consider this in the 
context of the warped decretion disc model, this delay can 
be explained as the time needed for the matter captured 


JD2400000+ 

X-ray flux 
(Crab units) 

<£> d 
102. 8 

42533 

2.8 

0.81 

42614 

0.2 

0.60 

42724 

0.3 

0.67 

42829 

0.3 

0.69 

43288 

0.5 

0.15 

43508 

0.5 

0.29 

43617 

0.2 

0.35 

43732 

0.7 

0.47 

43951 

0.2 

0.60 

44522 

1.5 

0.15 

44952 

0.2 

0.33 

45290 

0.2 

0.62 

45515 

0.5 

0.81 

45619 

0.8 

0.82 

45732 

0.8 

0.92 

46736 

0.8 

0.68 

47625 

0.6 

0.33 

49403 

1.4 

0.62 


Table 1 . Julian dates, X-ray fluxes and optical phases of 
the outbursts with F(2 — 10)keV > 0.2 Crab. X-ray data 
are from Giovannelli & Graziati (1992) and Finger et al. 
(1996). Giant outbursts of 1975 and 1980 are included. 


from the Be-star disc - disturbed during the close pas¬ 
sage of the neutron star - to travel to the vicinity of the 
NS companion, where it loses its angular momentum, falls 
onto the NS, and causes an X-ray flare. It can be noted 
that during 2001 two occasions of favourable X-ray and 
optical phases will be in June and October; unluckily, the 
former date corresponds to the seasonal gap of observa¬ 
tions. 


4. Conclusion 

An analysis of the uniform photometric data set obtained 
in the period 1983-1998 has allowed a confident separation 
of the periodic constituent in the light curve of the high- 
mass X-ray binary, A0535+26/V725 Tau. The parameters 
of this periodic component and its link with the phase 
of activity of the optical component allow us to suggest 
precession of an accretion disc around the neutron star or a 
warped equatorial disc around a Be star as the most likely 
mechanisms. At this point, both models seem to be viable; 
meanwhile, the analysis of already existing spectral data 
could be helpful, in the sense that if the warped disc model 
does reflect reality, one might expect to see V/R and EW 
variations corresponding to the processional motion. 

Within both models, we do not expect that substan¬ 
tial X-ray outbursts to occur during ascending parts of the 
large-scale optical light curve. Moreover, X-ray outbursts 
tend to occur at specific phases of the 102/8 optical light 
curve. Taken together, both effects can explain the “miss¬ 
ing outburst” phenomenon. 
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Our results lead us to suppose that in other similar 
systems we can hope to distinguish the minor photometric 
variations close to, but not necessarily coincident with, the 
binary rotation period. In future papers of this series, we 
plan to apply the techniques described here to the analysis 
of the light curves of similar systems, such as X Per and 
V635 Cas. 
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